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Abstract 

A realistic two-center shell model for fusion is proposed, which is based on two 
spherical Woods-Saxon potentials and the potential separable expansion method. 
This model describes the single-particle motion in a fusing system. A technique for 
calculating stationary diabatic states is suggested which makes use of the formal def- 
inition of those states, i.e., they minimize the radial nonadiabatic coupling between 
the adiabatic states. As an example the system ^^O -|- ^'^Ca ^^Ni is discussed. 
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1 Introduction 



The two-center shell model (TCSM) is a basic microscopic model to de- 
scribe the single-particle (sp) motion during a heavy ion collision near the 
Coulomb barrier. The sp energies and wave functions are obtained solving 
the Schrodinger equation with a phenomenological two-center mean field po- 
tential. This approach was first introduced by the Frankfurt school [1,2] and 
most of its applications so far in fission and fusion have been based on a double 
oscillator potential [3] . Improved versions of TCSM based on oscillator poten- 
tials have been developed for dealing with very asymmetric fission [4] and with 
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asymmetric fragmentations involving deformed fragments [5]. The two-center 
problem for realistic finite depth potentials has been solved with the wave 
function expansion method (usual diagonalization procedure) in Refs. [6,7]. 
The authors of Ref. [6] made use of two spherical Woods-Saxon potentials for 
describing the polarization of sp states in peripheral transfer reactions like 
^°Ca + ^^O and ^°^Pb + ^^C, while in Ref. [7] two deformed Gaussian poten- 
tials were applied for the fusion of a deformed and arbitrarily orientated ^^C 
projectile on the ^^O target. 

In the present work the two-center problem is solved using two Woods-Saxon 
(WS) potentials and the potential separable expansion (PSE) method pro- 
posed by Revai [8,9,10,11]. In contrast to Refs. [10, 11], where the PSE tech- 
nique was employed within the two-center framework to describe the sp motion 
in peripheral coUisions (the WS potential parameters were kept fixed) like ^^O 
_l_ i6q^ i7q _|_ i3q i7q _|_ i6q^ extend the potential ansatz for large 

overlap between the nuclei by conserving the volume of the system during the 
amalgamation process. For the sake of simplicity, we will use two spherical WS 
potentials, whose parameters in a realistic calculation have to be adjusted so 
as to fit the experimental sp spectra around the Fermi level for the separated 
nuclei and for the spherical compound nucleus. For intermediate internuclear 
distances these parameters can be interpolated as we will discuss below. The 
present TCSM has the following advantages compared to the traditional one 
based on two oscillator potentials [2]: (i) very asymmetric reactions can be 
described as those used to form superheavy elements (e.g., ^°Ca -|- ^°^Pb) 
without assumptions regarding the shape of the potential barrier between the 
fragments, (ii) the relative level positions in the colliding nuclei arc correct, 
(iii) the two-center potential barrier between the nuclei is realistic, (iv) the 
Coulomb interaction for protons can be explicitly included, (v) the sp wave 
functions have a correct asymptotic behaviour [10], and (vi) the continuum 
states like Gamow resonances can be included [11] which is relevant in reac- 
tions involving drip-line nuclei. It is worth mentioning that the PSE method 
has also been used in conjunction with both the one-center problem [12,13] 
and the two-center problem [14] for nuclear structure studies. 

In sect. 2, the PSE method to solve the adiabatic two-center problem is briefiy 
presented along with the method to ensure the volume conservation of the fus- 
ing system. The adiabatic or Born-Oppenheimer approximation is based on the 
idea that the relative motion of the nuclei is much slower than the sp motion, so 
at each fixed internuclear distance the molecular sp wave functions diagonalize 
the two-center sp hamiltonian. These states are called adiabatic sp states and 
the nucleons occupy the lowest sp energy levels in the two-center potential. 
The coupling between the adiabatic sp states, called nonadiabatic coupling, is 
induced by the kinetic energy operator related to the relative motion of the 
nuclei. When the adiabatic sp states are well separated, this coupling may be 
negligible and the adiabatic approximation is adequate. Whenever the energy 
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separation between two adiabatic sp states becomes small or vanishes, the 
coupling may be large and the adiabatic approximation breaks down. In this 
case, the sp state is a linear combination of adiabatic states. In such a situa- 
tion a useful approach, called stationary diabatic approximation, achieves the 
separation between the sp motion and the relative motion of the nuclei. The 
stationary diabatic sp states minimize the radm/ nonadiabatic coupling [15]. In 
the last part of sect. 2, we describe a technique to calculate the stationary dia- 
batic sp states following this formal definition. The convenience of the diabatic 
sp states in the entrance phase of heavy ion coUisions around the Coulomb 
barrier was discussed by Gassing and Norenberg in Ref. [16]. In contrast to the 
adiabatic sp motion, nucleons with a diabatic motion do not always occupy 
the lowest energy levels, but remain in their diabatic levels during a collective 
motion of the system. The use of the diabatic states accounts for the main part 
of the coherent coupling between collective and intrinsic degrees of freedom. 
Very recently in Ref. [17] these states were applied for the modeUing of the 
compound nucleus formation in the fusion of heavy nuclei. In sect. 3 numerical 
results for ^^O -|- ^°Ca ^^Ni are discussed, while conclusions are drawn in 
sect. 4. 



2 The TCSM 



We will follow the method proposed by Gareev et al. in Ref. [10] and also used 
by Milek and Reif in Ref. [11], which is based on the expansion of potentials in 
terms of, e.g., harmonic oscillator functions (PSE method), i.e, V ~ Xlij \i > 
Vij < j\ = Vsep- The accuracy of the results obtained depends only on the 
accuracy of the expansion of the WS potentials. The Schrodinger equation 
with approximate potentials Vgep is solved exactly. 

The finite depth nuclear potential belongig to each fragment (s = 1,2) is 
chosen to be the following spherical WS with a spin-orbit term 

W = -V../-M + lA,(-|-)Vo.lf(l.s), (1) 



where h/moc = 0.21 fm, / and fso are the same function, but with different 
parameters, i.e., fl,„){r) = {1 + exp[(r - ^o(so))/«o(5o)]}~^- ^or protons, the 
Coulomb potential Vq^i is taken to be that of a uniformly charged sphere with 
charge ZgC {Zg being the total charge of each fragment) and the radius R^, 
which is added to expression (1). The potentials (1) are placed at the position 
Rg in the center of mass system and thus the two-center potential reads as 

V = e-'^'^ Vi e^^i^ + V2 e^^^^ (2) 
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where k — hT^p is the sp wave-number operator. Each potential (1) is then 
represented approximately (Kjgp) within a truncated sp harmonic oscillator 
basis, {\v >,z/ = 1,...A^}, with the spin-angular part coupled to the total 
angular momentum j with projection m, e.g., in the momentum representation 

\v >= \nljm >= gni{k) ■ [r'Y,(k) ® x\ (3) 

The harmonic oscillator basis has the advantage that all matrix elements 
needed can be calculated analytically (e.g., see Appendix A in Ref. [11]) and 
this basis adopts essentially the same mathematical form in momentum and 
coordinate representations. In the coordinate representation (spherical coordi- 
nates) , the radial part of the harmonic oscillator wave function can be written 
as < r\nl >— h^^J"^ (l)ni{x) , where hose is the oscillator length and the adi- 
mensional x — h~}^ r. In the momentum representation, < k\nl >= gni{k) — 
(— 1)" /Iq/^ <PniiO being the adimcnsional ^ = hose k. The phase factor 
in the definition of the basis set (3) removes the imaginarity of gni{k). The 
function (p^i is the reduced oscillator function, i.e., 

where Ljj+^/^(,^^) is a Laguerre polynomial. According to Ref. [19], the oscillator 
lenght hose associated with each basis set can be calculated with the expression 
hl^^ = 0.84 • ■ ^;V6^ being As the mass number of each fragment. 

Volume conservation. To describe the amalgamation process of the two nuclei, 
the potential parameters have to be interpolated between their values for the 
separated nuclei and the spherical compound nucleus. The parameters can 
be correlated by assuming the condition that the volume enclosed by certain 
equipotential surface Vq of the two-center potential V is conserved for all values 
of i?, i.e., 

j dv e(Vo -V)^ vo, (5) 

where is the Heaviside step function, Vq the fixed equipotential surface cor- 
responding to the Fermi level of the fused system, and Vq the volume enclosed 
by the equipotential surface Vq. As the two-center potential V in (5), we use 
Vi + V2 where Vs refers to the first term of (1), i.e., the spin-orbit part is 
neglected. For the sake of simplicity, we will assume that the nuclear shape is 
the same for neutrons and protons. 

The parameters V — {Vos, -Ro(so)' "^0(50)' ^s, Rl, Zs} are interpolated by means 
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of the unknown function y{R) [7] 



V{K) = V{R = 0) [1 - y{R)] + V{R ^ oo) y{R), (6) 

where y{R) satisfies the boundary conditions y{R = 0) = and y{R oo) = 
1. The function y{R) G [0,1] is numericaUy calculated by inserting (6) into 
(5) and requiring the conservation of the volume. The parameters V{R 
oo) are unambiguously found fitting the experimental sp spectra around the 
Fermi level for the separated nuclei, but this is not the case for the individual 
fragment potentials at i? = where only the sp spectra of the compound 
system are known. At = we adopt that -Ro{so); "^0(50)' ^"^^ ^1 ^'^'^ the same 
as those for the compound nucleus, but Vqi + V02 = and Zi + Z2 = Zcn- 
A new parameter a — Vqi/V^^^ — Zi/Zcn is introduced and it determines 
the values of Vqs and at i? = 0. This a e (0, 1) parameter is chosen in such 
a way that the splitting of the sp levels for small R is similar to that in the 
Nilsson model. 

Adiabatic states. In the present method the two-center problem is solved in 
the momentum representation. The approximate potential Vf^p reads as 

Kep= E \sly>V:^<s^^\, (7) 



where z^, /i = {nljm}. The number of basis states included in the expansion 
(7) is defined by Imax (number of partial waves in which the potential acts) 
and Umax (the number of separable terms in each partial wave) . The values of 
Imax and Umax are determined by the convergence of the sp energies which is 
accelerated using the technique of the Lanczos cr-factors [18] . For bound states 
the formal solution of the Schrodinger equation is as follows 

\ip>=Go{E)V\ip>, G,{E)^{E-p^)-\ (8) 



where Go is the Green operator for the free sp motion. Inserting (2) and (7) 
into (8), the expression (8) reads as 

N 2 

\V >= Go{E) Y: E e-R-^^ \su >, (9) 



where 



A,^ =< s//|e^^^^|(^ > . (10) 
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The coefficients Agp, are the amphtude of the molecular sp wave function > 
with respect to the moving basis states located at R^. Multiplying (9) from 
the left by < sfile^^"^ and with (10), the following set of linear algebraic 
equations for the coefficients Ag^ is obtained: 



E Ei^ss'S,,^ - E < si,\Go{E)e'^^^-^^'^'\s'i^ > Ki,] = 0. (11) 



The solvability condition of the system (11) is that its determinant vanishes. 
This determines the energy eigenvalues E which appear as a parameter in the 
matrix elements containing Gq. With the eigenvalues E, the eigenvectors > 
are obtained solving the system (11) for the coefficients and requiring 
the normalization of the state vectors \(p >. These states are the adiabatic sp 
states. 

The quantum numbers denoting the molecular sp states arise from the sym- 
metry properties of the two-center problem. Since the potential (2) is invariant 
around the axis connecting the centers of both potentials (R = Ri — R2), the 
projection of the total angular momentum j along this axis denoted by Q is 
conserved. This axis is chosen as the quantization axis and the sp levels are 
labelled by Q and an additional index indx that distinguishes different states 
with the same Q. The index indx could refer either to the asymptotically good 
quantum numbers (nlj) of the state to which the level goes for large R or just 
to an index n counting the levels. In the intrinsic (or rotating) molecular ref- 
erence frame (Ri = R and R2 = 0) the matrix elements containing Go are as 



follows (j = (2j + 1)V2) 

s = s': <nljn\Go{E)\n'l'j'n' >=Sii,6jfSnn' < nl\Go{E)\n'l' >, (12) 

s^s' : 



N 



2 



N 



L 




<nl\Go{E) jL{kR)\n'l' >, 



(13) 



where 




(14) 
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<,.|G„(E)|,.T>^-^/<fM^,V*, 7^^-^. (15) 

^ 

/ nn (T7\ ■ (1 DM '7'^ J g*ni{k)gn'v{k)jL{kR) 2,, 

< nl\Go{E) jL{kR)\nl >= — — J ^ k dk. (16) 

^ 

In (16) Jl refers to the spherical Bessel function. The matrix elements (13) are 
obtained using (i) the expansion of e*-^''^''^ in terms of the spherical functions 

oo L 

giR,,,k^4^^ E ^''JL{kRss')YlM{^ss')yLMO^), (17) 

L=0 M=-L 

and (ii) general properties of the angular momentum algebra [20]. 

The sp states \indx Vl > calculated with (11)-(16) refer to adiabatic molecular 
sp states in the intrinsic (or rotating) reference frame, where the amplitude of 
the molecular wave functions solving (11) are denoted by Bgisn with /3 = {nlj}. 
Please note that these amphtudes depend only on the relative distance R. They 
are independent of the orientation angles R of the molecular axis, since these 
are no longer contained in the matrix elements (12)-(16). The rotation of the 
molecular axis with respect to a space-fixed axis is then included in (9) by 
means of the Wigner DI^q{R) rotation matrices, so the \indx fl > adiabatic 
molecular sp states in the laboratory reference frame read as follows 



(18) 



Stationary diabatic states. A natural way of introducing the stationary dia- 
batic sp states is to impose the condition that the radial coupling operator 
{—ihd/dR + A) be zero or negligibly small in these states [15], where A iden- 
tifies and cancels that spurious portion of the coupling induced by —ihd/dR 
that represents merely the translation of basis functions along with moving 
nuclei. Fortunately only inspection is needed to remove from the —ihd/dR 
matrix that part which corresponds simply to displacement, as we will show 
next. Then what remains after such background correction is the actual nona- 
diabatic coupling which vanishes for large internuclear distances. The diabatic 
representation is obtained by diagonalization of this coupling in a truncated 
set of the adiabatic sp states (18) with the same quantum number Q, i.e., 
< indx fl\{— ihd/dR + A)\indx' fl >. Only couplings between neighbouring 
levels {indx' — indx ± 1) are included. In defining the matrix of the radial 
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nonadiabatic couplings that is diagonalized, we neglect the couplings whose 
strength is smaller than a threshold value '-fthr- The value of 'jthr used in the 
calculations should lead to the isolation of the strongest radial nonadiabatic 
couplings. Strong radial nonadiabatic couplings are expected to occur at the 
avoided crossing of the adiabatic molecular sp levels. 

The adiabatic molecular sp states (18) can be re- written as 

\indx n >= Go(£^md.,n){E Cf.e"^^' 11^^ > + E ^2.1'^^ >>' (19) 

u=l v=l 



where the coefficients = ^up'n^sp'n- Here, we identify u = {/3m} and 
the Wigner rotation matrices are now included in the basis {siy >. The radial 
derivative d/dR of (19) is as follows 



d 



iy=l u=l 

^1 f)r<^ ^2 



+ Go(£;.„..,n){E ^e-I^^ \lu > + E ^ |2- >} 
+ Go{E,r,a.,n) E ^i^^^ > ■ (20) 



The last term of (20) is responsible for the unphysical effect (i.e. non van- 
ishing asymptotic couplings owing to the displacement of the moving basis) 
on the radial nonadiabatic transitions discussed above and will therefore be 
removed. Consequently, from expression (20) one can observe that the ra- 
dial nonadiabatic coupling vanishes for i? — > 00 because dGo{Eindx,n)/dR = 
-dEindx,n/dR ■ Gl{Eindx,n) and dC^^/dR when the molecular sp 
states approach the fixed sp states of the separated nuclei. In Appendix A we 
give expressions for the calculation of the actual radial nonadiabatic coupling 
matrix elements including the first four terms of (20) only. 

After diagonalizing the actual radial nonadiabatic coupling we obtain the uni- 
tary transformation matrix Sfj between the adiabatic and diabatic represen- 
tations, where i,j refer to the values of indx. The diabatic sp wave functions 
and energy levels are obtained as follows 

\i ^ >diab= E^i]- ^ >adiab, (21) 
j 

rpdiab \ ^ i „f! |2 rpadiab /'oo^ 

^ia = 2^ . {22) 

j 
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For a given Q the adiabatic sp levels E^q'^'' as a function of R (correlation 
diagram [22]) show avoided crossings because of the noncrossing rule by von 
Neumann and Wigner [21], while the diabatic levels Ef^^ reveal real cross- 
ings. The diabatic states (21)-(22) can cross each other because they are not 
solutions of an eigenvalue problem. 

There arc approximating methods [22] for calculating diabatic states, namely 
the maximum symmetry method and the maximum overlap method. These 
techniques were applied in Ref. [23] to obtain diabatic states with a TCSM 
based on harmonic oscillators. In Refs. [17,24] we made use of the maximum 
symmetry method to calculate diabatic potential energy surfaces of near sym- 
metric reactions. Other methods have been developed in quantum chemistry 
for the calculation of such a diabatic basis, e.g., see Ref. [25]. 



3 Numerical results 

As an example of the realistic TCSM discussed above we present calculations 
for neutrons and protons in the reaction ^^O -|- "^''Ca ^^Ni. The potential 
parameters for the separated nuclei and the compound nucleus as well as 
the maximal number of nodes n^nax and partial waves l^ax included in the 
harmonic oscillator basis for converged values of the sp energies are given in 
Table 1. These parameters approximately reproduce the experimental sp levels 
around the Fermi energy, which was shown in Refs. [6,26]. 

3. 1 Neutrons 

Fig. 1 shows the function y{R) derived from the condition of volume con- 
servation for different values of the a parameter which determines both the 
partition of the compound nucleus potential depth into two parts belonging 
to each fragment and their charges. The numerically calculated function y{R) 
is used in expression (6) to interpolate the two-center potential parameters. 
It is observed that (i) y{R) decreases from 1 around the contact configuration 
{R ~ 7.3 fm) up to for the spherical compound nucleus (i? = fm), and (ii) 
y{R) very strongly depends on the a parameter at small radii R, while this is 
not the case around the contact point. The strong effect of the a parameter 
on the neutron sp spectra at small distances R can be seen in Fig. 2. 

In Fig. 2 we show normalized neutron sp levels with Q = 1/2 as a function 
of the quadrupole deformation parameter j3 of the ^^Ni compound nucleus for 
two values of the a parameter, namely a = 0.4 (right panel) and a = 0.8 
(middle panel). Moreover, these level diagrams are compared to the levels of 
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Table 1 

Potential parameters for the separated nuclei and the compound nucleus: The po- 
tential depth is given in MeV, the radii and diffusenesses in fm, while the strength 
A of the spin-orbit potential is adimensional. The radii -Ro(so) = ^o(so)^^''^ being A 
the mass number of the nucleus. See text for further details. 



Neutrons 



Nucleus 


^0 




ao 


A 




a-so 


Rc 




Imax 


16Qa 


74.90 


0.97 


0.60 


19.50 


1.10 


0.55 




4 


3 


^OCa'^ 


51.40 


1.28 


0.76 


24 


1.37 


0.55 




4 


3 




46.20 


1.31 


0.62 


18.73 


1.31 


0.62 




4 


3 










Protons 










16 QO 


79.10 


0.95 


0.59 


19.15 


1.10 


0.55 


2.77 


4 


3 


40Ca" 


70.80 


1.05 


0.67 


19.22 


1.05 


0.55 


3.83 


4 


3 


56]\^j6 


53.70 


1.24 


0.63 


13.97 


1.31 


0.63 


4.74 


4 


3 



Prom Ref. [6], ''Prom Ref. [26] 



the Nilsson model [27] (left panel). The normalization is given by the first 
Nilsson-level and the f3 value is obtained from the quadrupole moment of the 
nuclear shape. From this figure it is concluded that a realistic value of a is 
about 0.8 because, as expected, the behaviour of the two-center sp levels is 
similar to that in the Nilsson model. In the following we will use a — 0.8. 

Fig. 3 shows the central part of the neutron two-center potential along the 
z-a.xis for fixed values of the distance R between the fragments. The nuclear 
shape associated with each potential is also included in this figure. Near the 
contact configuration {R — 7 fm) the interfragment potential barrier practi- 
cally disappears and, therefore, most of the nucleons may move in the whole 
volume of the system. From distances i? < 4 fm the two-center potential essen- 
tially adopts a one-center character. It is worth mentioning that the two-center 
potential barrier is naturally formed in the present model by the superposition 
of the tail of the two WS potentials. 

Fig. 4 shows the adiabatic (left panel) and diabatic (right panel) neutron 
levels with = 1/2 as a function of the internuclear distance R. For small 
and large radii we can observe the shell structure of the compound nucleus 
and the separated fragments, respectively. For distances 3 fm < i? < 7 fm sp 
levels are strongly polarized and the shell structure of the separated fragments 
is practically dissolved. Here sp excitations and transfer processes essentially 
occur around the avoided crossing of the molecular adiabatic sp levels due to 
the effect of the radial nonadiabatic couplings on the sp motion. Most of the 
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1.2 




R(fm) 



Fig. 1. The function y{R) derived from the condition of volume conservation for 
different values of the a parameter. This is used to interpolate the neutron and 
proton two-center potential parameters as a function of the radius ill between the 
nuclei. See text for further details. 

avoided crossings in the adiabatic level diagram turn into real crossings in the 
diabatic level diagram. This happens where the radial nonadiabatic coupling 
between the adiabatic sp levels shows huge peaks which can be seen in Fig. 5. 

Fig. 5 shows the absolute value of the radial nonadiabatic couplings between 
some neighbouring adiabatic levels of Fig. 4 (left panel) as a function of the 
internuclear radius: levels 3-4 (solid curve), levels 5-6 (dotted curve) and levels 
8-9 (dashed curve). For large distances the couplings vanish as expected. The 
more bound the levels are, the smaller is the radius R in which the coupling 
vanishes, e.g., comparing levels 3-4 (solid curve) and levels 8-9 (dashed curve). 
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'"Ni (Q = 1/2) 




0.000 0.025 0.050 0.025 0.050 0.025 0.050 

P 

Fig. 2. Normalized neutron levels with = 1/2 as a function of the quadrupole 
deformation P of the compound nucleus '''^Ni for different values of the a parameter: 
a = 0.4 (right panel) and a = 0.8 (middle panel). The left panel shows the levels 
of the Nilsson model. See text for further details. 

The huge peaks localize the position of an avoided crossing where a nucleon 
transition (Landau-Zener effect) may occur with a large probability. These 
peaks are the key to calculating the diabatic states which are obtained by 
minimizing those strong couplings only. The threshold value for the strength 
of the couplings was jthr = 1-4. It is important to stress that strong transitions 
between molecular sp levels as that between the levels 8-9 (dashed curve) 
at i? ~ 8 fm can also occur for distances larger than that in the contact 
configuration. In these situations Landau-Zener nucleon transitions may be 
reflected in the excitation function of direct reaction processes [3] . 
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T 




z(fm) 

Fig. 3. The central part of the neutron two-center potential along the z-axis for 
fixed radii R between the nuclei. The nuclear shape associated with each potential 
is also included. See text for further details. 

In Fig. 6, the adiabatic neutron levels correlation diagram including all bound 
states with different fl values is presented, i.e., f2 = 1/2 (solid curves), fl — 3/2 
(dashed curves), Q = 5/2 (dotted curves) and Q = 7 /2 (dashed-dotted-dotted 
curve). The states with the same Q value arc coupled by the radial coupling 
discussed above, while those states with fl values differing by one unit in fl 
(e.g., fl = 1/2 and fl = 3/2) are coupled in the laboratory reference frame 
by the so-called rotational (or Coriolis) coupling which is maximal at the real 
crossing of these states [3] . 

Both the adiabatic and the diabatic sp states with different fl values, as men- 
tioned above, are coupled in the laboratory system by the rotational coupling. 
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I ' I ' I ' I ' I ' I ' I ' ' I ' I ' I ' I ' I ' I ' I ' 

2 4 6 8 10 12 2 4 6 8 10 12 

R(fm) 



Fig. 4. Adiabatic (left panel) and diabatic (right panel) neutron levels with f2 = 1/2 
as a function of the radius R between the nuclei. Some adiabatic levels are labeled 
by 3-9. See text for further details. 

The rotational coupling (~ j-I, where j is the nucleon total angular momentum 
and I is the total angular momentum of the system) appears in a non central 
collision due to the rotation of the molecular internuclear axis with respect 
to a space-fixed axis. Please note that the sp wave-functions ((18) and (21)) 
depend on the Wigner rotation matrices. One can expect that the rotational 
coupling (which increases with decreasing internuclear radii R, i.e., ~ 
is important in light ion collisions owing to the small moment of inertia of 
the system, but it is not the case in a heavy ion reaction. In fusion of light 
systems (where the fusion process is already determined by the penetration 
of the external Coulomb barrier, i.e., small overlap between the nuclei) the 
relevant radial region for the rotational coupling is governed by the usually 
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Fig. 5. Absolute value of the radial nonadiabatic coupling between some neighbour- 
ing adiabatic levels of Fig. 4 (left panel) as a function of the radius R between the 
nuclei. See text for further details. 



strongly absorptive nucleus-nucleus potential. Terlecki et al. [28] have found in 
the reaction ^'^C+^'^C that the effect of the rotational coupling on the inelastic 
excitation cross sections is much smaller than the effect of the radial coupling. 
We would like to add that Imanishi and von Oertzen have introduced in Ref. 
[29] the so-called Rotating Molecular Orbitals (RMO) in the standard cou- 
pled reaction channels formalism in order to avoid numerical problems with 
the rotational coupling at small internuclear radii R. To our knowledge, there 
has been no systematic study of the effect of rotational couplings on fusion 
and scattering cross sections. Further work on this issue is required. 
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Fig. 6. 

Fig. 7. The adiabatic neutron levels correlation diagram including all bound states 
with different O values. See text for further details. 

3.2 Protons 



Since we have assumed that the nuclear shape is the same for neutrons and 
protons, the same function y{R) of Fig. 1 (a = 0.8, dashed-dotted curve) is 
apphed in interpolating the proton two-center potential parameters. In Fig. 
7 like in Fig. 3, the central part of the proton two-center potential along the 
2;-axis is shown for different fixed internuclear distances R. Apart from (i) the 
positive tails of the two-center potential for large z values due to the Coulomb 
interaction and (ii) a higher barrier between the fragments at i? = 7 fm, 
the change of the two-center potential with R is similar to that in Fig. 3 for 
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z(fm) 

Fig. 8. The same as Fig. 3, but for protons. See text for further details. 



neutrons. Fig. 8 shows the adiabatic (left panel) and diabatic (right panel) 
correlation diagram of the proton levels with Q = 1/2. The same features as 
those observed for neutrons in Fig. 4 can be seen here. In order to isolate 
the strong peaks of the proton radial nonadiabatic coupling, a threshold value 
Ithr = 1-7 for the strength of the couplings was needed in this case to calculate 
the diabatic levels. Fig. 9 is like Fig. 6, but for protons. The level diagram 
shows similar features as those in Fig. 6, only the levels are shifted up due to 
the Coulomb interaction. 
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Fig. 9. The same as Fig. 4, but for protons. See text for further details. 
4 Concluding renicirks 

A realistic TCSM for fusion has been proposed, which is based on two spher- 
ical WS potentials and the PSE method. This model describes the sp motion 
in a fusing system. The adiabatic and diabatic molecular sp states for neu- 
trons and protons in the reaction ^^O -|- ^°Ca — > ^^Ni have been calculated. 
A technique to compute the stationary diabatic states has been introduced, 
which is based on the minimization of the strong radial nonadiabatic coupling 
between the adiabatic states. The only difference between the stationary dia- 
batic sp basis and the adiabatic one is that the diabatic basis minimizes the 
strong radial nonadiabatic coupling localized at some of the avoided crossings 
between the adiabatic levels. Whether it is more convenient to use either the 
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Fig. 10. The same as Fig. 6, but for protons. See text for further details. 

adiabatic or the diabatic basis essentially depends on the bombarding energy 
[16]. It is important to add that the dynamical radial nonadiabatic coupling is 
proportional to the ion-ion radial velocity R which modulates the stationary 
(or static) couplings calculated in the present work. 

At very low incident energies (e.g., energies below the Coulomb barrier) one 
can expect that the real sp motion is well described by the adiabatic repre- 
sentation. The dynamics of the sp motion becomes even simpler when the dy- 
namical nonadiabatic couplings may be neglected and so the nucleons remain 
at the lowest energy levels. When the incident energy increases, the dynamical 
couplings arc more and more strengthened. In this case, nuclcon transitions be- 
tween the adiabatic molecular sp levels occur. The adiabatic basis along with 
all nonadiabatic couplings should then be used for understanding the real sp 
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motion. If the incident energy reaches a certain vahic (e.g., 1 MeV/u above 
the Coulomb barrier for heavy symmetric systems [16]) at which the dynami- 
cal nonadiabatic coupling is very strong (probability of the nucleon transition 
is close to one), then the diabatic sp basis approximates the real sp motion. 
At this energy regime, the use of the diabatic basis is certainly very useful to 
simplify the dynamical calculation. It is worth mentioning that the diabatic sp 
motion is expected to be realistic in the entrance phase of heavy-ion reactions 
[16]. The adiabatic and diabatic representations can be used as sp bases in 
a coupled channels calculation within a molecular quantum-mechanical for- 
mulation of heavy ion collisions [3,29]. Works regarding the apphcation of 
the present TCSM in understanding the formation mechanism of superheavy 
elements are in progress. 
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A Actual radial nonadiabatic coupling matrix elements 



The matrix elements < indx fl\d/dR\indx' Q > can be calculated as follows 
including only the first four terms of (20) (for clarity Q will be removed, and 
new a and /3 indexes will refer to indx and indx', respectively) 



d 

< a\——\(3 >— terml + term2 + ■ — h termS, 
oR 



(A.l) 



where 



Ni Ni 

terml = Yl < 1^' 

u'=l u=l 

N2 N2 

term2 = £ ^ C^^, < 2u' 

v'=l u=l 

Ni Ni 

terms =J2Y1 ^'iV < 1^' 

u'=l u=l 

N2 N2 

termA = £ E ^1' < 2^' 

I/' = l u=l 

Ni N2 

termb = J2 J2 ^i^' < 1^' 

u'=l i/=l 



OR ' 

dC2u 
OR ' 



GoiE^f-^^ \2u > C,., 



G^{E^)Go{Ep)e^^ \2v > 



(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 



20 



termG = E E < '^A Go{E^)Go{Ep)e-'^^ \\v > (A.7) 

„'=1 u=l 

term! = E E ^iV < 1^1 ^°(^") Qif ^ '^^ > (^"S) 



i/'=i i/=i 

7V2 iVi 



terms = E E ^1' < 2^1 go(^a) ^ > ^1- (A-9) 



=1 i/=i 



The matrix elements contained in ierml andierm2 (i.e., < su'l G'o(£'a)G'o (£'/?) > 
with s — 1,2) are obtained using both the relation 

Go{Ea)Go{Ep) = —^-—[Go{E^) - Go{Efs)], E^ ^ Ep (A.IO) 



and the matrix elements given by expression (12). 

The matrix elements of termS and termA are calculated applying the relation 

< siy'l go(^a) g^ "^ Isiy >= < Go{E^)Go{E^) Isiy > .(A.ll) 

The rest of the matrix elements involved in termb, ■ ■ ■ , termS are obtained 
making use of (i) the relation (A.IO), (ii) the matrix elements calculated with 
expression (13), and (iii) the derivative rule used in the relation (A.ll). 
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